Extended surfaces modulate and can catalyze hydrophobic effects 
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Abstract 

Interfaces are a most common motif in complex systems. To understand how the presence of 
interfaces affect hydrophobic phenomena, we use molecular simulations and theory to study hydra- 
tion of solutes at interfaces. The solutes range in size from sub- nanometer to a few nanometers. The 
interfaces are self-assembled monolayers with a range of chemistries, from hydrophilic to hydropho- 
bic. We show that the driving force for assembly in the vicinity of a hydrophobic surface is weaker 
than that in bulk water, and decreases with increasing temperature, in contrast to that in the bulk. 
We explain these distinct features in terms of an interplay between interfacial fluctuations and ex- 
cluded volume effects — the physics encoded in Lum-Chandler- Weeks theory [J. Phys. Chem. B 103 
4570-4577 (1999)]. Our results suggest a catalytic role for hydrophobic interfaces in the unfolding 
of proteins, for example, in the interior of chaperonins and in amyloid formation. 
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Hydrophobic effects are ubiquitous and often the most significant forces of self-assembly 
and stability of nanoscale structures in liquid matter, from phenomena as simple as micelle 
formation to those as complex as protein folding and aggregation fD, S| . These effects depend 
importantly on lengthscale pHS]. Water molecules near small hydrophobic solutes do not 
sacrifice hydrogen bonds, but have fewer ways in which to form them, leading to a large 
negative entropy of hydration. In contrast, hydrogen bonds are broken in the hydration of 
large solutes, resulting in an enthalpic penalty. The crossover from one regime to the other 
occurs at around 1 nm, and marks a change in the scaling of the solvation free energy from 
being linear with solute volume to being linear with exposed surface area. In bulk water, 
this crossover provides a framework for understanding the assembly of small species into a 
large aggregate. 

Typical biological systems contain a high density of interfaces including those of mem- 
branes and proteins, spanning the entire spectrum from hydrophilic to hydrophobic. While 
water near hydrophilic surfaces is bulk-like in many respects, water near hydrophobic sur- 
faces is different, akin to that near a liquid- vapor interface [5H5]. Here, we consider how 
these interfaces alter hydrophobic effects. Specifically, to shed light on the thermodynamics 
of hydration at, binding to, and assembly at interfaces, we study solutes with a range of sizes 
at various self-assembled monolayer interfaces over a range of temperatures using molecular 
simulations and theory. 

Our principal results are that although the hydration thermodynamics of hydrophobic 
solutes at hydrophilic surfaces is similar to that in bulk, changing from entropic to enthalpic 
with increasing solute size, it is enthalpic for solutes of all lengthscales near hydrophobic 
surfaces. Further, the driving force for hydrophobically driven assembly in the vicinity of 
hydrophobic surfaces is weaker than that in bulk, and decreases with increasing temperature, 
in contrast to that in bulk. These results suggest that hydrophobic surfaces will bind to and 
catalyze unfolding of proteins, which we predict are relevant in the formation of amyloids 
and the function of chaperonins. 

A. Models 

Molecular simulations: We simulate the solid-water interfaces of self-assembled mono- 
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FIG. 1. Size dependent hydrophobic hydration at interfaces, (a) A schematic of a cuboidal cavity 
(green) at the SAM- water interface. The SAM head-groups (black and white), alkane tails (grey), 
and water (red and white, partially cut out for clarity) are shown, (b) A typical configuration of 
the model membrane, color-coded by its distance from the model surface (grey), (c) Important 
volumes in estimating the free energy ^'^^ of emptying the probe volume V (green) using the 
theoretical model. The region above the membrane is the volume B (blue), and the intersection 
of V and B is v (dark green), (d) Lengthscale dependence of the cavity hydration free energy per 
unit area, fi^^/A, in bulk water and at interfaces, at T = 300 K, obtained from MD simulations, (e) 
fi^^/A, estimated using the theoretical model, near surfaces with different attraction strengths ij. 
(f) Connecting the microscopic binding free energy of a cavity to an interface, to the macroscopic 
surface wettability. The cos 9 values were obtained from MD simulations of a water droplet on 
SAM surfaces ([9]). Lines are predictions using Eq. ^ with size dependent 7lv taken from (d). 

layers (SAMs) of surfactants [Fig. [l]^a)] with a range of head-group chemistries, from hy- 
drophobic (-CH3) to hydrophilic (-0H) [71 19]. To study the size dependence of hydration 
at interfaces, we selected cuboid shaped [L x L x W) cavities, with thickness W = 0.3 nm, 
and side L, varying from small values comparable to the size of a water molecule to as large 
as ten times that size. Thicker volumes will show qualitatively similar behavior, but will 
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gradually sample the "bulk" region away from the interface. 

Theoretical model: To rationalize the simulation results and obtain additional physical 
insights, we developed a model based on Lum- Chandler- Weeks (LCW) theory [4J. LCW 
theory incorporates the interplay between the small lengthscale gaussian density fluctua- 
tions and the physics of interface formation relevant at larger lengthscales, and captures the 
lengthscale dependence of hydrophobic hydration in bulk water. Near hydrophobic surfaces, 
it predicts the existence of a soft liquid-vapor-like interface, which has been confirmed by 
simulations [BHH]- 

We model this liquid-vapor-like interface near a hydrophobic surface, as an elastic mem- 
brane [Fig. [l|^b)], whose energetics are governed by its interfacial tension and the attractive 
interactions with the surface. The free energy of cavity hydration, /i''^, is related to the 
probability of spontaneously emptying out a cavity shaped volume, V. Such emptying can 
be conceptualized as a two-step process in which interfacial fluctuations of the membrane 
can empty out a large fraction of V in the first step, with the remaining volume v emptied 
out via a density fluctuation [Fig. [l]^c)]. When v is small, the probability that it contains 
waters is well-approximated by a Gaussian [71 [101 E] • The cost of emptying v can then be 
obtained from the average and the variance of number of waters in v, which are evaluated by 
assuming that the water density responds linearly to the surface-water adhesive interactions. 

We tune the strength of the model surface-water attraction, f/(r), using a parameter rj, 
where t] ^ 1 corresponds to the hydrophobic -CH3 SAM-like surface, with higher values 
representing increasingly hydrophilic surfaces. The representation of hydrophilic surfaces 
in our theoretical model lacks the specific details of hydrogen bonding interactions (e.g., 
between the hydrophilic -OH SAM surface and water), so comparisons between high-?] model 
surfaces and hydrophilic SAM surfaces in simulations are qualitative in nature. Equations 
that put the above model on a quantitative footing are given in the Appendix and the details 
of its exact implementation are included as Supplementary Information. 



B. Size-dependent hydrophobic hydration at, and binding to interfaces 

Fig.[T](d) shows the excess free energy, /i®^, to solvate a cuboidal cavity at temperature T = 
300 K, divided by it's surface area {A = + ALW). i^^/A can be thought of as an 
effective surface tension of the cavity-water interface. In bulk water, this value shows a 
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FIG. 2. Temperature dependence of in bulk water and at SAM-water interfaces for (a-b) large 
{L = 3.0 nm) and (c-d) small {L = 0.5 nm, L = 0.75 nm) cavities, obtained from simulations (a,c) 
and from the model (b,d) of Eq. ([s]). 

gradual crossover with increasing L, as expected |ll[l2]- Fig-[l]^d) also shows the lengthscale 
dependence of fi^^/A for solvating cavities in interfacial environments. Near the hydrophilic 
OH-terminated SAM, the behavior is similar to that in bulk water. However, with increasing 
hydrophobicity of the interface, the size dependence of fi'^^/A becomes less pronounced and is 
essentially absent near the -CH3 surface, suggesting that hydration at hydrophobic surfaces 
is governed by interfacial physics at all lengthscales. 

Fig. [l](e) shows the analogous solvation free energies predicted using the theoretical model. 
The essential features of solvation next to the SAM surfaces are captured well by this 
model. This is particularly true for the hydrophobic surfaces (with t] around 1), where the 
potential U{r) closely mimics the effect of the real SAM on the adjacent water, and the 
agreement between theory and simulation is nearly quantitative. For the more hydrophilic 
SAMs, the comparison is qualitative, because the simple form for U (r) does not represent 
dipolar interactions well. 

Fig. [l|d) also indicates that n'^^ becomes favorable (smaller) with increasing surface hy- 
drophobicity. The difference in fi^^ at an interface and in the bulk, Afi"^^ = fi'?^^ — /i^uiki 
quantifies the hydration contribution to the experimentally measurable free energy of bind- 
ing of solutes to interfaces. Because the solvation of large solutes is governed by the 
physics of interface formation, both in bulk and at the SAM surfaces, we can approxi- 
mate Afi"^^ = Ac(7sv — 7sL — 7lv), where Ac = L"^ is the cross-sectional area, 7 is the surface 
tension, and subscripts SV, SL, and LV, indicate solid-vapor, solid-liquid, and liquid-vapor 
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interfaces, respectively. Using Young's equation, 7sv = 7sl + 7lvcos^, we rewrite 



A/i''" = -Ac7Lv(l-cos^), (1) 

where 9 is the water droplet contact angle on the solid surface. 

Although Eq. ([T]) is strictly valid only for macroscopic cavities, it can be applied to 
sufficiently large microscopic cavitities with a lengthscale-dependent surface tension, 7lv 
(~ ^'^^Yk/A). Indeed, lines in Fig. [l|^f) predicted using Eq. ([T]) are in excellent agreement 
with simulation data, and indicate that the strength of binding increases with surface hy- 
drophobicity as well as with solute size. These results establish a connection between the 
microscopic solute binding free energies to interfaces and the macroscopic wetting properties 
of those interfaces. This connection provides an approach to characterize the hydrophobicity 
of topographically and chemically complex interfaces, such as those of proteins [121 E] . 

C. Temperature dependence of hydration at interfaces 

The differences between the cavity hydration at interfaces and in bulk are highlighted 
most clearly in the T-dependence of /i'^^, which characterizes the entropic and enthalpic con- 
tributions to the free energy. For small solutes in bulk, the entropy of hydration is known 
to be large and negative [121 US] , which reflects the reduced configurational space available 
to the surrounding water molecules. In contrast, for large solutes the entropy of hydration 
is expected to be positive, consistent with the temperature dependence of the liquid-vapor 
surface tension [TT] . Fig. [2]^a) shows that jj,^^ of large cuboidal cavities (L = 3 nm) in bulk 
water indeed decreases with increasing temperature, although the corresponding hydration 
entropy per unit surface area (25 J/mol/K/nm^) is lower than that expected from the tem- 
perature derivative of surface tension of water (about 90 J/mol/K/nm^ [l^J). We note that 
solvation entropies in SPC/E water obtained using NPT ensemble MD simulations are known 
to be smaller than experimental values by about 20% [H]. Additionally, the cavity- water 
surface tension and its temperature derivative for these nanoscopic cavities are expected to 
be smaller than the corresponding macroscopic values [8] . 

Fig. [2|^a) also shows that for large cuboidal cavities (L = 3nm), jjL^^ decreases with 
increasing temperature not only in bulk water and near the hydrophilic (-OH) surface, but 
also near the hydrophobic (-CH3) surface, indicating a positive entropy of cavity formation. 
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Thus, in all three systems, the thermodynamics of hydration of large cavities is governed 
by interfacial physics. Although the values of fi^^{L = 3nm) at 300 K are rather large 
(569kJ/mol in bulk water, 565kJ/mol at the -OH interface, and 167kJ/mol at the -CH3 
interface), their variation with temperature shown in Fig. |2|a) is similar in bulk and at 
interfaces. 

Fig. shows that this same phenomenology is captured nearly quantitatively by the 
theoretical model. In the model, the cavity hydration free energies have large but ather- 
mal contributions from the attractions between water and the model surface. The main 
temperature-dependent contribution to fi^^ is the cost to deform the liquid-vapor-like inter- 
face near the surface to accommodate the large cavity. Since the necessary deformation is 
similar, regardless of the hydrophobicity of the surface, the variation of /i*^^ with temperature 
is similar as well. 

Fig. |2|^c) shows the temperature dependence of fi^^ for small cavities (L = 0.5 nm) in bulk 
and at SAM-water surfaces. In bulk water, /i*^^ increases with temperature, and yields an 
entropy of hydration of roughly — 25 J/mol/K, characteristic of small lengthscale hydropho- 
bic hydration. This negative value is consistent with those calculated for spherical solutes 
of a similar volume p^J. With increasing hydrophobicity, the slope of the fi'^^ vs T curve 
decreases and becomes negative, indicating a positive entropy of cavity formation near suf- 
ficiently hydrophobic surfaces. Near the most hydrophobic surface (-CH3), the entropy of 
hydration of this small cavity is +30 J/mol/K. 

Fig. |2](d) shows that the same phenomenon is recovered by the theoretical model, though 
the correspondence is clearest at a slightly larger cavity size (L = 0.75 nm). Near hydrophilic 
model surfaces, the interface is pulled close to the surface by a strong attraction, so it is 
costly to deform it. As a result, the cavity is emptied through bulk-like spontaneous density 
fluctuations that result in a negative entropy of hydration of small cavities. In contrast, near 
a hydrophobic surface, the interface is easy to deform, which provides an additional mecha- 
nism for creating cavities. In fact, this mechanism dominates near sufficiently hydrophobic 
surfaces, and since the surface tension of water decreases with increasing temperature, so 
does 11^^. Hence, even small cavities have a positive entropy of hydration near hydrophobic 
surfaces. The continuous spectrum of negative to positive solvation entropies observed in 
Figs. |2](c-d) is thus revealed to be a direct consequence of the balance between bulk-like 
water density fluctuations and hquid-vapor-like interfacial fluctuations. 
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FIG. 3. Lengthscale dependence of the excess solvation entropy per unit surface area for (a) 
cavities in bulk water and at the -CH3 and -OH SAM-water interfaces, and (b) cavities in the 
model of Eq. ([s]) near surfaces of different attraction strengths, rj. 

Fig. |3|^a) shows that near the hydrophobic CHs-terminated SAM, cavity hydration en- 
tropies per unit area, S'^'^/A, are positive and essentially constant (about 30 J/mol/K/nm^) 
over a broad range of cavity sizes. In contrast, in bulk water, S^^/A depends on L, and 
changes from large negative to positive values with increasing L. The lengthscale at which 
entropy crosses zero, Lg, can serve as a thermodynamic crossover length. In bulk water, 
Lg ^ 1.8 ± 0.2 nm. The behavior of S'^^/A is qualitatively similar at the -OH surface, with 
Lg ^ 1.3 ± 0.4 nm. Although the numerical value of Lg may depend on the shape of the 
cavity and on solute-water attractions for non-idealized hydrophobes, the trend in entropy 
should not. 

Fig. shows that our implementation of LCW ideas recovers many of the ob- 

served trends, with solvation entropy being everywhere positive for the smallest attraction 
strength rj, and a thermodynamic crossover length of just under Inm emerging for the 
more hydrophilic model surfaces, similar to that in bulk water. Nevertheless, the agreement 
between Figs. [3|a) and (b) is somewhat qualitative, mostly as a result of the crude form of 
U (r) used to model hydrophilic surfaces. 



D. Thermodynamics of binding to, and assembly at hydrophobic surfaces 

In the preceding sections, we have examined the hydration behavior of single, isolated, 
idealized cavities near flat surfaces and in the bulk. We now consider the consequences of our 
observations on hydrophobically driven binding and assembly, summarized schematically in 
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Fig-i 

Fig. |4] indicates that while the binding of both small and large solutes (or aggregates) 
to hydrophobic surfaces is highly favorable, their thermodynamic signatures are different. 
Binding of small solutes is entropic and becomes more favorable with increasing temperature, 
whereas binding of large solutes is enthalpic and depends only weakly on temperature. 

Fig. |4] also highlights the differences in the thermodynamics of hydrophobically driven 
assembly at interfaces and in bulk, inferred from our lengthscale dependence studies. In 
bulk, the solvation of many small, isolated hydrophobes scales as their excluded volume. 
Accommodating small species inside the existing hydrogen-bonding network of water im- 
poses an entropic cost, so the solvation free energy increases with increasing temperature. 
When several small hydrophobes come together, water instead hydrates the aggregate by 
surrounding it with a liquid-vapor-like interface. The corresponding solvation free energy 
scales as the surface area and decreases with increasing temperature. 

Thus, the driving force for assembly of n small solutes (each of surface area Ai, volume Vi 
and solvation free energy of /^i^uik) i^^^ ^ large aggregate (with surface area and volume 
nvi) in bulk water is well-approximated by 

A/^buik = 7buik^„ - n/i^^buik = [Cri-'/' - l]^/ii"buik, (2) 

where C ~ (Tbuik^^i''^/ A^i^uik) ^^'^ 7buik is a curvature-corrected effective surface tension [top 
curve of Fig. [l]^d)]. As the surface tension decreases with increasing temperature, so does the 
free energy to hydrate nanometer-sized aggregates. However, the free energy to individually 
hydrate the small solutes increases with temperature, resulting in a larger driving force for 
assembly. Conversely, while the driving force for assembly, A/i^^jj^, is large and negative 
(favorable) at ambient conditions, it decreases in magnitude with decreasing temperature 
[upper portion of Fig. |4] , and can even change sign at a sufficiently low temperature. When 
adapted to particular systems, Eq. ^ can, with remarkable accuracy, explain complex 
solvation phenomena like the temperature-dependent aggregation behavior of micelles [12] 
and the cold denaturation of proteins [2]. 

In the presence of a hydrophobic surface, on the other hand, we have found that interfacial 
physics dominates at all lengthscales [Fig. |2]^a-d) and Fig. [3]^a-b)]. As a result, the driving 
force for assembly at interfaces, Ayufj^^, does not scale as in Eq. ([2]), but is instead given by 

= 7int(A„ - riA,) ~ [n-'/' - l]n/x^,^i,„ (3) 
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where 7int is the effective surface tension at the interface, [the lower curves of Figs. [T]^d- 
e)]. Since 7int decreases with increasing temperature [Figs. |2](a-d)], so does the hydration 
contribution to the driving force for assembly at a hydrophobic surface, in contrast to that 
in bulk. 

^ Assembly in bull< 




FIG. 4. Schematic illustrating the thermodynamics of binding and assembly. The points represent 
free energies of solvating small objects individually (left) and in the assembled state (right), in 
bulk (top) and at a hydrophobic interface (bottom), at a lower (blue, Tl) and a higher (red, Th) 
temperature near ambient conditions. Assembly: The driving force for assembly at hydrophobic 
interfaces is smaller than that in bulk. It is enthalpic and decreases with increasing temperature, 
unlike in bulk. Binding: The driving force for binding small objects to a hydrophobic surface 
increases with temperature, so it is entropic. For large objects, it is enthalpic. 

The free energy barrier between disperse and assembled states is also expected to be 
very different in bulk and near hydrophobic surfaces. In bulk, the dispersed state has no 
liquid-vapor-like interface whereas the assembled state does. The transition state consists of 
a critical nucleus of hydrophobic particles that nucleates the liquid-vapor-like interface. The 
nucleation barrier can be high, and dominates the kinetics of hydrophobic collapse of ideal- 
ized hydrophobic polymers [2DH22] and plates [2Bj. In contrast, we expect aggregation near 
hydrophobic surfaces to be nearly barrierless, since an existing liquid- vapor-like interface is 
deformed continuously between the disperse and assembled states. 

Finally, and most importantly, we find that for large aggregates, the driving force of 
assembly is weaker near interfaces than in bulk. In the limit of large n, the terms n//^^ 
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dominate both at interfaces and in bulk (Eqns. ^ and (|3|), and the results in Fig. [T|^d) 
show that fxl^-^, < /i^^^uik- 

The nontrivial behavior of the driving forces and barriers to assembly at interfaces should 
be relevant in biological systems where hydrophobicity plays an important role. Experiments 
have shown that hydrophobic surfaces bind and facilitate the unfolding of proteins, including 
those that form amyloids |2^] - [26] . Our results shed light on these phenomena and suggest 
that large hydrophobic surfaces may generically serve as catalysts for unfolding proteins 
[27] . via solvent-mediated interactions. Indeed, simulations show that the binding of model 
hydrophobic polymers to hydrophobic surfaces is accompanied by a conformational rear- 
rangement from globular to pancake-like structures [23 . Such conformations can further 
assemble into secondary structures, such as /3-sheets [2SH2ZI 121], and we predict that the 
solvent contribution to this assembly at the hydrophobic surface will be governed by inter- 
facial physics. This implies that manipulating the liquid-vapor surface tension, either by 
changing the temperature or by adding salts or co-solutes, will allow one to manipulate the 
driving force for assembly. 

We further speculate that the catalysis of unfolding by hydrophobic surfaces may play a 
role in chaperonin function [30j • The interior walls of chaperonins in the open conformation 
are hydrophobic and can bind misfolded proteins, whereupon their unfolding is catalyzed [311 
[32] • Subsequent ATP-driven conformational changes render the chaperonin walls hydrophilic 
[301 [31] . As a result, the unfolded protein is released from the wall, as the free energy for a 
hydrophobe to bind to a hydrophilic surface is much lower than that to bind to a hydrophobic 
one [Fig.[gd)]. 

Our results also provide insights into the interactions between biomolecules and nonbi- 
ological hydrophobic surfaces, such as those of graphite and of certain metals, which have 
been shown to bind and unfold proteins [3Sl [31] . Such interactions are of interest in diverse 
applications including nano-toxicology [35j and biofouling [34J. 

Collectively, our findings highlight that assembly near hydrophobic surfaces is different 
from assembly in bulk and near hydrophilic surfaces. Experimental measurements of the 
thermodynamics of protein folding have been performed primarily in bulk water [36]. Al- 
though many experiments have probed how interfaces affect protein folding, structure and 
function [211 131], to the best of our knowledge, there are no temperature-dependent thermo- 
dynamic measurements of self-assembly at interfaces. We hope that our results will motivate 
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such measurements. 



APPENDIX 

Simulation details: Our simulation setup and force fields are similar to that described 
in Refs. [3, |9]. Simulations were performed in the NVT ensemble with a periodic box 
(7nmx7nmx9nm) that has a buffering liquid- vapor interface at the top of the box, for 
reasons explained in Ref. [8J. It has been shown that free energies obtained in the above 
ensemble are indistinguishable from those obtained in the NPT ensemble at a pressure of 
Ibar [38j. We have chosen the SPC/E model of water [39] since it adequately captures 
experimentally known features of water, such as surface tension, compressibility, and local 
tetrahedral order, that play important roles in the hydrophobic effect [5j. Electrostatic 
interactions were calculated using the particle mesh Ewald method |10j, and bonds in water 
were constrained using SHAKE [H]. Solvation free energies were calculated using test 
particle insertions |12] for smaller cavities [L < Inm), and the indirect umbrella sampling 
(INDUS) method [HI EH] for larger cavities. 

Theoretical Model: We model the liquid-vapor-like interface near hydrophobic surfaces 
as a periodic elastic membrane, z = h{x,y), with an associated Hamiltonian, H[h{x,y)]: 



H[h{x,y)] 



x,y 



l\Vh{x,y)\'+ / p,f/(r) 

, ^ J z=h{x,y) 



(4) 



Here, 7 is the experimental liquid-vapor surface tension of water, pe is the bulk water density, 
and U{r) is the interaction potential between the model surface and a water molecule at 
position r = (x, y, z). The square-gradient term in Eq. (|4]) accurately captures the energetics 
of interfacial capillary waves only for wavelengths larger than atomic dimensions (Fig. |5]), 
so we restrict h{x,y) to contain modes with wavevectors below 27r/9A. At any instant in 
time, part of V can be empty due to an interfacial fluctuation. The number of waters in the 
remaining volume, v, fluctuates, and we denote by Pv{N) the probability that v contains 
waters. We thus estimate the free energy for emptying V completely to be 

p<'-(V) = -fceTln j mZ-ie-''^[''("'^)lp,(0), (5) 

where Z = J Vh exp{— /3i^[/i(x, y)]} is the partition function of the membrane. The vol- 
ume V depends on the interfacial configuration h{x,y), i.e., v = v[h{x,y)]. 
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It is known that Pv{N) is well-approximated by a Gaussian when v is small [TJ [101 
[TT] . If water were far from liquid- vapor coexistence, then Pv{N) would also be close to 
Gaussian for arbitrarily large v. The fact that water at ambient conditions is near liquid- 
vapor coexistence, and that there is a liquid-vapor-like interface near the SAM, is captured 
by the additional interfacial energy factor exp{—(3H[h{x, y)]} in Eq. ([s]). The net result 
is that the thermal average of Eq. ([s]) is dominated by interface configurations where v is 
small, so that even at ambient conditions, we can approximate 

where (N), is the average number of waters in v and a, = {{SN)'^), is the variance. We 
estimate these by noting that the solvent density responds linearly to the attractive po- 
tential, f/(r), in the volume occupied by the water, B, depicted in Fig. ^c) [101 HSl HI]. 
Hence, 

{N),^pev- [ [ x(r,r')/3t/(r0, 

JrGv J r'eB 

o'v ^ / x(r, r'), where 
X(r,r') = pi6{r - r') + p|[^(|r - r |) - 1]. 

Here, g{r) is the oxygen-oxygen radial distribution function of water |3H]. 

The surface-water interaction is modeled by a potential, U{r), that closely mimics the 
attractive potential exterted by the -CH3 SAM on water: 

U{r) = t/waii(r) + r/?7head(r) + t/taii(r). 

The first term, f/waii(r), is a sharply repulsive potential in the region z < Rq that captures 
the hard-core exclusion of a plane of head groups at z = with hard-sphere radius Rq. The 
second term, f/headli"), captures the head group-water interaction, modeled as a plane of 
OPLS/UA CH3 Lennard- Jones (LJ) interaction sites [49] of area density //head at 2; = 0, and 
is scaled by 77. The final term, ?7taii(i"), similarly captures the alkane tail-water interaction, 
modeled as a uniform half-space of OPLS /UA CH2 L J interaction sites of volume density ptaii 
at a distance ( below the head groups. The parameters Rq, (, phead and ptaii are dictated 
by the geometry of the SAM (See SI for details). 
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FIG. 5. Power spectrum of the instantaneous liquid- vapor interface at T = 300 K. A liquid- vapor 
interface was simulated using a 24 x 24 x 3nm^ slab of SPC/E water in a periodic box of size 
24 X 24 X 9nm^ and the instantaneous interface configuration, h{x,y), and its Fourier transform, 
/i(k), were evaluated as in Ref. |l6]. The power spectrum of our simulated instantaneous interface 
is good agreement with the capillary-wave theory prediction ((|/i(k)p) ~ l//3^k'^) for wavevectors 
smaller than ~ 27r/9 A. For larger wavevectors, the power spectrum is sensitive to molecular detail, 
i.e., the coarse-graining length ^ used to define the intrinsic interface, as expected |47]. Fitting 
the ^ = 2.0 A data in the range 0.01 A ^ < k < 0.3 A ^ yields 7 = 62.0 it 0.5 mJ/m^, in reasonable 
agreement with the experimental value of 72mJ/m^ and some simulated values of the SPC/E 
surface tension {e.g., 63.6 ± 1.5mJ/m^ [S])) but not others {e.g., 52.9mJ/m^ [17]). 
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SUPPLEMENTARY INFORMATION 



Here, we describe the details of the solvation model used in the main text. 



1. Interface description 

We describe the liquid-vapor-like interface next to the model surface by a periodic height 
function /i(a), with a = {x,y) and —D/2 < x,y < D/2. This function is sampled discretely 
at a resolution A, at points satisfying 

/ . .X D D 

a={n^A,nyA), -^<n^,ny<—. 

This results in N"^ discrete sampling points {a}, with N = D/ A. In the following, sums 
over a denote sums over these A^^ sampling points. We have used D = 60A and A = 1 A. 
The discrete variables {hg,} represent the interface height at each sample point a, so that 

ha = /i(a), for a = {nxA,nyA). 

This notation clearly distinguishes between the N"^ height variables and the continuous 
height function /i(a) that they represent. 

The discrete Fourier transform of {/ia} is denoted by {/ik}, and is defined at wavevec- 
tors k = (2,71/ L) {nix, my), with —N/2 < nix, my < N/2. We use the symmetric normaliza- 
tion convention throughout for Fourier transforms. 



2. Energetics 

The essential property of the liquid- vapor-like interface is its surface tension, which results 
in the following capillary-wave Hamiltonian [50] for a free interface, 

a k 

where Vh^ is a finite-difference approximation to V/i(a). 

Using an appropriate definition of an instantaneous water- vapor interface ^6], the power 
spectrum of capillary waves in SPC/E water has been found to agree with the spectrum 
predicted by the above Hamiltonian for wavevectors smaller than about 27T/i, but is sub- 
stantially lower for higher wavevectors (Fig. 6 of main text). This result is consistent with 
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the liquid- vapor-like interfaces being sensitive to molecular detail at high wavevectors |17] . 
At T = 300 we have found that £ ^ 9 A. We thus constrain all Fourier components /ik to 
be zero for high k, i.e. 

K = 0, |k| > 27i/£. (6) 

In our model, the liquid-vapor-like interface interacts with a model surface via a potential 
that depends on {ha}. As discussed below, it is also convenient to introduce additional 
umbrella potentials to aid in sampling. The Hamiltonian of the interface subject to this 
additional potential energy H'[{ha}] is 

k 

When expressed as a function of the Fourier components {/ik}, we denote the Hamiltonian 
by -?/[{/ik}] and the external potential by i7'[{/ik}], so that 



3. Dynamics 

We calculate thermal averages of interface configurations by introducing a fictitious 
Langevin dynamics and replacing thermal averages by trajectory averages. We first assign 
a mass per unit area fi to the interface. The Lagrangian in real space is 



^[{/^a,M] = VE^a-^[{M]- 



The corresponding Lagrangian in Fourier space is 

k 

Since all are real, the amplitudes of modes k and — k are related, /ik = h*_^^. Taking this 
constraint and Equation (|6| into account, the Euler-Lagrange equations yield equations of 
motion in Fourier space. To thermostat each mode, we add Langevin damping and noise 
terms. The final equation of motion has the form 

fiA^ = -jA'lkfh^ - ^^JjiM _ r^n^ + (|k| < 2n/i), (8) 
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The Langevin damping constant r] is chosen to decorrelate momenta over a timescale r, so 
1] = yuA^/r. The zero-mean Gaussian noise terms {^k(^)} have variance such that 

mt)Ut')) = 2r^k^TSit-t'). 

As with /ik, satisfy the related constraint .^k = 'C*k- Hence, for k = 0, the noise is 
purely real and its variance is twice that of the real and imaginary components of all other 
modes [51]. 

We propagate these equations of motion using the Velocity Verlet algorithm. At each 
force evaluation, we use a Fast Fourier Transform (FFT) to calculate {ha} from {/ik}- We 
then calculate dH'[{ha}]/dha in real space and perform an inverse FFT to obtain the force 
dH'[{hi^}]/dhi^ on mode hy^ due to H'[{ha}]- We then add the forces due to surface tension, 
Langevin damping and thermal noise, as in Eq. (|8]). 

For the Velocity Verlet algorithm to be stable, we choose a timestep equal to 1/20**^ of the 
typical timescale of the highest-frequency mode of the free interface. At = 2^\//i^V7- To 
equilibrate the system quickly but still permit natural oscillations, we choose the Langevin 
damping timescale so that r = lOOAt. Finally, we choose a value of close to the mass of 
a single water layer, /i = lOOamu/nm^. 

This interface dynamics is entirely fictitious. However, it correctly samples configurations 
of the interface Boltzmann-weighted by the Hamiltonian i^[{/ia}]- This is true regardless 
of the exact values of fi, At and r, so our choices have no effect on the results in the main 
text. We have simply chosen reasonable values that do not lead to large discretization errors 
when solving the system's equations of motion. 

4. Surface-interface interactions 

The liquid- vapor-like interface interacts with the model surface via a potential H!,^^f[{ha}]- 
In the atomistic simulations, the SAM sets up an interaction potential U (r) felt by the atoms 
in the water molecules. Below, we use the notation r and (x, y, z) interchangeably. To model 
this interaction potential, we smear out the atomistic detail of the SAM and replace it with 
three elements: 

• A uniform area density //head of Lennard- Jones sites (with length and energy scales 
(Thead and ehead) ill the z = Q plane to represent the SAM head groups. 
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• A uniform volume density ptaii of Lennard- Jones sites (with length and energy scales 
(Ttaii and etaii) in the half-space z < to represent the SAM tail groups. 

• Coarse-graining the head-group atoms into a uniform area density results in a softer 
repulsive potential allowing the interface to penetrate far deeper into the model surface 
than would be possible in the actual SAM. To rectify this, we apply a strongly repulsive 
linear potential in the half-space z < Rq, where Rq is the radius of the head group's 
hard core. The repulsive potential is chosen to be 1 ksT when 1 nm^ of interface 
penetrates the region z < Rohy a "skin depth" 5. 

The head groups are thus modeled by the following potential acting on a water molecule at 
position r: 

/oo /"OO 
dx' / dy' liLj (|r - r'|; ehead, o-head) |^,=o' 
-oo J — oo 

where ■ULj(r;e,(T) = 4e[((T/r)^^ — (o"/r)®] is the Lennard-Jones pair potential. Similarly, the 

effect of the tail groups is captured by 

/oo /'OO 
dx' / dy' / dz'-ULj (|r - r'|;etaii,o-taii) . 
oo J —oo J —oo 

Finally, the repulsive wall is modeled by the potential 

U^an{x,y,z < Ro) = 2p~^ ■ (1 /cBT/lnm^)(i?o - ^)/^, 

where = 0.03333 A ^ is the number density of liquid water. 

These smeared interaction potentials depend only on not on x or y. As described in the 
main text, we also scale the head- group interaction by a parameter 77. Putting everything 
together, we obtain an explicit expression for the surface-interface interaction potential, 

^surf[{M]=P^A2E^surf(^a), 
a 

where 

/•oo 

^surf(^a) = / d^r7C/head(^) + t^tail(^) + ^wal^^^), 
*/ /la 

^Trehead/^headC^head [^(^^head/^)^ " | (c^head/^)^] 
= ' +7retailPtailO-t\il [^((^head/[^ + C])^ " \{.CrheaJ[z + C])^] , Z > Rq, 

KnriiRo) + • (1 kBT/l nm2)([i?o - z\/5)\ z < i?o. 

To model the -CH3 SAM in this paper, we chose the following values for the parameters 
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The head groups are modeled as OPLS united-atom CH3 groups interacting with 
SPC/E water, so ahead = 3.5355 A and etaii = 0.68976 kJ/mol. 



• The tail groups are modeled as OPLS united-atom CH2 groups (sp^-hybridized) inter- 
acting with SPC/E water, so ataii = 3.5355 A and ctaii = 0.5664 kJ/mol. 

• The tail region is inset from the plane of the head groups by a distance equal to a CH2- 
CH3 bond length (1.50 A), minus the van der Waals radius of a CH2 group (I.9525A), 
so C = -0.4525 A. 

• The head group density is known from the atomistic SAM geometry to be /^head = 
0.0462 A^^. The mass density of the SAM tails was estimated to be 935kg/m^ [7], 
resulting in a CH2 group number density of ptaii = 

0.0402 A-3. 

• The equivalent hard sphere radius Rq of a -CH3 group at room temperature was 
estimated to be 3.37 A [Sj. It has a small temperature dependence, which we neglect. 

• The wall skin depth 6 was set to 0.1 A, which is small enough so that the repulsive 
potential is essentially a hard wall ai z = Rq, but large enough that we can propagate 
the interfacial dynamics with a reasonable timestep. 

5. Umbrella sampling 

Calculating fi'^^{V) from Equation (2) of the main text as a thermal average {Py{0)) over 
Boltzmann-weighted configurations of /i(a) is impractical for large V. The configurations 
that dominate this average simply have a vanishingly small Boltzmann weight. To solve 
this problem, and in analogy to what we do in atomistic simulations, we perform umbrella 
sampling on the size of the sub-volume v of the probe cavity V that is above the interface. 

We begin by defining the volume V corresponding to a probe cavity of dimensions L x 
L X W as the set of points satisfying \y\ < L/2 and Rq < z < Rq + W. We then 
define f [{/la}] as the size of the sub- volume of V that is above the interface. Using umbrella 
sampling and the multistate Bennet acceptance ratio method (MBAR) [52], we calculate the 
probability distribution for v, P{v), down to f = 0. To do this, we use quadratic umbrellas 
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defined by a center v and width {5vY, which result in the addition to the Hamiltonian of 

During each umbrella run, we also record the configurations {/ia} which yield each observed 
value of V. We then approximate the right-hand side of Equation (2) in the main text by 
summing over these configurations with appropriate weights, and obtain 

where, as in the main text, the term Pv{0) depends on the interface configuration {ha}, and 
the sum {/ia} is over all interface configurations in all the different umbrellas. To evaluate 
Pv{0), we implement discrete versions of the integrals defining {N)y and as was done in 
Ref. 
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